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ABSTRACT 



We use different assumptions for the combination of 
dust extinction and star formation up to redshifts z 
~> 5, as well as detailed modeling of supernova prop- 
erties, to estimate the number of supernovae that 
should be observable with various instruments, in 
particular the NGST. In the model we use realistic 
light curves and spectral shapes that evolve with time 
for the different types supernovae. 

We find that the NGST should be able to detect sev- 
eral tens of core collapse SNe in a single detection if 
the observational limit in the range 1-5 fim is ~ 1 
nJy and the field is 16 arcmin 2 . We also estimate the 
observable number of Type la supernovae. Due to 
the time delay between the formation of the progeni- 
tor star and the explosion of the supernova, there are 
additional aspects that have to be considered when 
interpreting observational rates of these supernovae. 

Key words: high redshift supernovae. 



1. INTRODUCTION 



The evolution of the star formation rate is reflected in 
the cosmic supernova rate (SNR). It should therefore, 
in principle, be possible to use SN observations to 
distinguish between star formation scenarios. Even 
more important, core collapse SNe, i.e. Types II and 
Ib/c, provide a direct probe of the metallicity produc- 
tion with cosmic epoch. Due to their importance as 
probes for early cosmological epochs, studies of high 
z SNe is one of the main goals for the Next Genera- 
tion Space Telescope (NGST) (Stockman & Mather 
1996[), 



In this paper we present estimates for the expected 
number of observable supernovae and discuss some 
of the complications entering the ana lysis. Previous 
studie s include e.g. Madau et al. (1998 ), Madau 
(ll997j ) Ruiz -Lapuente & Canal ( |l998| t"\ 7 ungelson & 
Livio ( 1998| ). With respect to these our work dif- 



fers in that we include information about the light 
curve, which allows us to predict the simultaneously 
observable number of SNe. We divide the SNe into 
different types with maximum absolute magnitudes 



accordingly. We also calculate the counts for different 
filters and include information about the estimated 
redshifts of the detected SNe. 

Section 2 describes our model. Results are presented 
in section 3. In section 4 we discuss the influence of 
dust. Further discussion follows in section 5. Conclu- 
sions are given in section 6. We assume a cosmology 

with H = 50 km s _1 Mpc 3 and qa=0.5. 



2. THE MODEL 



A problem when using luminosity densities from high 
z galaxies, to calculate the SFR, is that these can 
be matched to wide range of SFRs, spanning from 
strongly peaked to flat or even increasing at z<^l, by 
adjust ing the assumed extinction due to dust. Madau 
(1997) finds, by using luminosity densities in three 
bands (UV, optical and NIR) that a universal extinc- 
tion Eb-v =0.1 with SMC-type dust makes the best 
fit between a proposed SFR and the observations. 
This SFR shows a pronounced peak at 1< z <2, see 
Figure |l|. This shape is predicted in a scenari o whe re 



galaxies form hierarchically (e.g. Cole at al. 1994). 



In this paper we calculate the supernova rates 
by using the star formation rate and extinction 
E B -v=0.1, derived by Madau ( |1997| ). We also com- 
ment on how a higher extinction affects the rates, 
and if it is possible to use supernova counts to make 
predictions about the amount of dust. 

Core collapse supernovae, which are thought to origi- 
nate from short lived massive stars, have an evolution 
that closely follows the shape of the SFR. Assuming 
an immediate conversion of these stars to supernovae 
renders a multiplicative factor between the star for- 
mation and the supernova rates. 



SNR =k x SFR, 



(1) 



where k depends on the IMF and the mass range of 
the progenitors. 

The uncertainties concerning the life-time of Type 
la progenitors make the relation between the rate of 
these supernovae and the star formation more diffi- 
cult to model. We treat the life-time r of the progen- 
itor WD as a free parameter which is allowed in the 
range 0.3< r <3 Gyr. 
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Figure 1. Supernova rates for the model with dust ex- 
tinction Eb-v =0.1. The solid line represent the rate of 
core collapse supernovae. Dashed, dotted and dash-dotted 
lines represent the rate of Type la supernovae with time 
delays t=0.3, 1.0 and 3.0 Gyr respectively. The rates of 
Type la SNe have been normalized so that the ratio be- 
tween the core collapse and the Type la rates matches the 
locally observed value 3.5. Solid line together with scales 
on the right y-axis show the star formation rates. 



In order to derive the number of SNe to different 
limiting magnitudes, the apparent magnitude of the 
SNe are calculated from 

m f {z,t,i) = Mf(t) + n(z) + Kf(z,t) + A g + < A t f > 

(2) 

Here Mt{t) is the absolute magnitude in a filter / 
at time t relative to the peak of the light curve, 
fi(z) is the distance modulus, Kf(z,t) gives the K- 
correction, A g is the Galactic absorption, and < 
Aij > is the absorption in a parent galaxy with in- 
clination i. The core collapse SNe are divided into 
five different types, i.e. Ib/c, IIP, IIL, SN1987A-type 
and Iln, with absolute magnitudes accordingly. We 
model the spectra of the different SN types by black- 
body curves with cutoffs at short wavelengths. The 
spectra are also set the evolve with time. 

An important feature of our model is that the SNe 
are distributed over the light curve, and are given ab- 
solute magnitudes accordingly. With this procedure 
we obtain the simultaneously observable number of 
supernovae. In order to actually detect the SNe, a 
second observation has to be made after an appro- 
priate time has passed. The effect of the time delay 
between the formation of the progenitor star and ex- 
plosion of Type la SNe shifts the peak of the rate of 
these SNe towards lower z than compared to the core 
collapse SNe. In Figure we show the rates for core 
collapse and Type la SNe. Three different values for 
the life-time of Type la progenitors are shown, r=0.3, 
1 and 3 Gyr. The rates of Type la's are normalized 
so that the fraction between core collapse and Type 
la rates matches the observed value at z=0. 



3. RESULTS 



3.1. Core Collapse Supernovae 



Figure g shows the predicted number of simultane- 
ously observable core collapse SNe per arcmin 2 in 
the M' (centered on 4.2/im with AA = ±0.8/im), 
K' (2.1/im), J (1.2/xm), I and R filters for different 
limiting magnitudes (given in the AB system). Also 
shown are the mean redshifts of the SNe at the dif- 
ferent limits and the mean redshifts in Am=0.5 mag- 
nitude bins for the K' and M' filters. 

The number of observable high z SNe increases con- 
siderably more at faint limiting magnitudes for the 
filters that probe longer wave-lengths. This K- 
correction effect is enhanced by fact that a majority 
of the core collapse SNe has a spectral cutoff at short 
wave-lengths, causing these SNe to drop out of the 
R, I, J, K' and M' filters at z ~0.9, -1.4, -2.7, -5.5 
and —10, respectively. 

The NGST is planned to have a detection limit of —1 
nJy in the J and K' bands and —3 nJy in the M' band 
for a 10 4 s exposure with S/N=10 and R=3 (NGST 
Simulation Group, 1998). These fluxes correspond to 
AB-magnitudes of 31.4 and 30.2, respectively. With 
these limits we find that the K' band yields the high- 
est predicted number of observable SNe. In a 4x4 ar- 
cmin field we predict — 47 simultaneously detectable 
core collapse SNe in the K' band. The numbers for 
the J and M' bands are — 28 and — 41, respectively. 
The mean redshift of the SNe are < z >— 1.9 in 
the K' band and < z >- 2.1 in the M' band. The 
mean redshift in the faintest observable magnitude 
bin (Am=0.5 mag) for both these filters are < z >— 
2.3. The number of simultaneously observable SNe 
with z > 2 is - 17 for both the K' and the M' fil- 
ter. At even higher z, the M' filter is preferred due 
to the effects of K-correction. For SNe with redshifts 
z > 5 and z>9we predict 1.1 and 0.09 SNe in the 
M' band, while the numbers for the K' band are 0.3 
and 0.02 (using field and detection limits as above). 
These high z estimates are based on a supernova rate 
derived by extrapolating the rate from z< 5. This 
makes the predicted numbers somewhat uncertain. 
However, the M' band should detect a factor > 3-5 
times more SNe than the K' band at these z, indepen- 
dent of the actual rate. We stress that it is necessary 
to use filters that sample wave-lengths long ward of 
the (redshifted) spectral cutoff of the SNe spectra in 
order to detect SNe at high z. Note that the number 
of observable SNe with z above — 2-4 may be in- 
creased by the ef fects of gravitational magnification 
(Marri & Ferrara |1998|). 



An HDF-like observation, covering 5 arcmin 2 with 
observational limits of I<>29 and R<^29, should con- 
tain — 1.2 and — 0.6 core collapse SN in each 
field. The ESO Very Large Telescope (VLT) reaches 
R-2 6.8 using FORS in 1800s exposure (Moorwood 
1994|) . With a field of 6.8' x 6.8', one should detect - 



.5 core collapse SNe. The mean redshift of these SNe 
is predicted to be < z >- 0.8. About 30% and 0.8% 
of these have z >1 and z >2, respectively. Reaching 
a magnitude deeper would raise the detected number 
to — 3 and increase the mean redshift to < z >~ 1. 
With this limit 40 % and 2 % of the SNe have red- 
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Limiting magnitude Limiting magnitude 



Figure 2. Left panel: the number of core collapse supernovae per arcmin 2 that can be detected down to different limiting 
magnitudes. Results for M , K , J, I and R bands are shown as full, dashed, dotted, dash-dotted and dash-triple- dotted 
lines, respectively. Right panel: the mean redshift of the SNe (same representation as left panel). Filled and open dots 
represents the redshifts in magnitude bins of Am— 0.5 for the M and K' bands. Note that AB-magnitudes are used. 



shifts z >1 and z >2. Using the VLT infrared ISAAC 
instrument, which is specified to reach K' >24.8 in 
a 2.5'x2.5' field (3600s exposure) (Moorwood |1994| ), 
will yield around 0.07 detected SNe per observation 
(< z >^ 0.8). Reaching a magnitude deeper in- 
creases the predicted number of observable SNe to 
0.4 per field (< z >~ 1). In the last example about 
47% and 4% of the SNe have z >1 and z >2, re- 
spectively. These estimates show that ground based 
instruments can not compete with the NGST when 
it comes to detecting high z supernovae. 

So far we have discussed the number of SNe that 
are simultaneously observable during one search (re- 
quiring two observations to actually detect the SNe). 
Prolonging the time coverage of a specific field in- 
creases the observed number of SNe. With 4-6 ob- 
servations per year, yielding a full coverage of the 
light curves, we estimate a total of ~ 70 SNe to be 
detected per field with the NGST in the K' filter. 

Compared to earlier published estimates we predict 
greater numbers for the count s. With the same ex- 
tinction, Madau et al. ( 1998|) predict ~7 SNe per 
NGST field per year in the range 2< z <4. Our cal- 
culations result in ^23 (both K' and M'). The main 
reason is that we include SNe observed over the whole 
light curve. 



3.2. Type la Supernovae 

Figure || shows the detectable number of Type la SNe 
down to different limiting magnitudes for r=0.3, 1 
and 3 Gyr in the different filters. Also shown is the 
mean redshift for the model with r=l Gyr. 



introduces a fairly large dispersion in the predicted 
number of observable SNe. For example, with the 
NGST K' detection limit ~ 12 SNe are predicted for 
the two low t values and ~ 7 for t=3 Gyr. This 
may seem as a possible observational probe to distin- 
guish between different progenitor models. However, 
in next section we show that this use is hampered by 
the uncertainties introduced when considering alter- 
native extinction models. 



4. DEPENDENCE ON DUST 



There are several claims that the extinction 
Eb-v=0.1, as argued for by Madau (1997), and used 
here, is too low. Observations of high z galaxies, 



as well as the f ar-infr ared extra-galactic backg round 



(e.g. Meurer 
Sa wicki & Yee 
al, 



1922 



1995 



Rowan- Robinson, 
Ellingson ct al 



et 



al 

19961 Soi 



1997 

fcr 



1997| ), indicate t 
■^0.3 could be more likely 



et 

at an ex- 
With 



1998| ; Burigana et al 
tinction of Eb-v 
a high extinction it is possible that the observed 
peak in the UV-luminosities from the high z galax- 
ies is illusive and that the star formation history 
is more compatible with a monolithic collapse sce- 
nario. To mim ic the case of a monolithic collapse, 
Madau ( 1997[ ) presents a test case where the ex- 
tinction is adjusted so that a proposed monolithic 
SFR matches the observations. In this model he 
finds an extinction law that increases with redshift 
as £B_ y =0.011(l+z) 2 - 2 . 

We have run our model with the monolithic extinc- 
tion law in order to examine how a higher absorption 
a high redshifts affects the predicted rates. 



These results show that a change in the delay time, r, 
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Limiting magnitude Limiting magnitude 

Figure 3. Left panel: the number of predicted Type la SNe per arcmin 2 in the M 1 , K 1 , J, I and R filters for different 
limiting magnitudes. The full, dashed and dotted lines represents T=0.3, 1.0 and 3.0 Gyr. The left hand scaling is valid 
for the M 1 counts. The counts in other filters are offset by factors 10 , 10~ 2 , 10~ 3 and 10~ A , respectively. Left panel: 
the lines show the mean redshifts of the SNe at the different limits for the T=l Gyr model. Filled and open dots show the 
mean redshift in magnitude bins (Am=0.5) for the Af and K bands, respectively. Note that AB-magnitudes are used. 



4.1. Core collapse SNe 

The rate of exploding high z core collapse SNe in- 
creases in proportion to the assumed increase in star 
formation in the high extinction model. Due to the 
effects of the dust, the observable number of SNe 
will, however, be similar to the numbers found in 
the low extinction model. This is hardly surprising, 
because the observed supernova rate is, except for 
the K-correction, proportional to the observed star 
formation rate, which (by assumption) is the same in 
the two models. With the NGST limits, the number 
of detected SNe in a field increases from ^47 in the 
low dust scenario to ~54 in the high dust scenario, 
using the K' band. 

Since the K-corrections increase with redshift, it may, 
however, be possible to use a high redshift subsample 
of the observed SNe in the K' or M' band as a probe 
for the amount of dust. Using NGST limits, and 
only counting SNe with z > 2 in the low and high 
dust models, results in ~ 17 and ~ 32 detections, 
respectively (in both bands). 



4.2. Type la SNe 

The increased star formation at high redshifts in the 
monolithic model makes the rate of high z Type la 
SNe sensitive to the choice of formation redshift, zf- 
Allowing the formation redshift to be in to range 5 
< zf < 9, together with the uncertainty in progenitor 
lifetime r, leads to a wide dispersion in the estimated 
counts of Type la SNe. Using the NGST limits re- 
sults in ^5-36 simultaneously observable Type la SNe 
in the K' filter. A better understanding of the prop- 



erties of the dust seems necessary if the rates are to 
be used as probes for different progenitor scenarios. 



DISCUSSION 



The actual amount of extinction has a large influence 
on the cosmic star formation rate above redshifts of 
the peak at z ~1.5. This affects the intrinsic number 
of supernovae and therefore also the production and 
enrichment of metals by core collapse SNe through- 
out history. The number of core collapse supernovae 
as observed from Earth is, however, less dependent 
on the amount of extinction, due to the fact that 
high intrinsic star formation rates are combined with 
high absorption as to match the same observational 
data. A more pronounced change in the derived rates 
would emerge if the observed UV luminosities from 
the high z galaxies under-predicts the star formation 
due to other effects than extinction. This could be 
the case if a large part of the star formation takes 
place in low surface brightness galaxies, or if selec- 
tion effects biases against the detectio n of g alaxies 
with increasing redshift (e.g. Hu et al. 
son|!998D. 



1998; Fergu- 



The rates of simultaneously detectable SNe results 
from events over the whole light curve. Using Type 
la's as standard candles requires observations at the 
peak of the light curve, i.e that a first detection is 
made at the rising part of the light curve. Our calcu- 
lations predict that the NGST will be able to simul- 
taneously detect ~ 0.6 Type la's on the rising part 
of the light curve, of which ~ 0.4 and ~ 0.1 has red- 
shifts z >1 and z > 2, respectively. We have used 
the hierarchical model with extinction £^s_y=0.1, 
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and progenitor life-time r=l Gyr, in these examples. 
The dependence of these numbers on the progenitor 
life-time is, however, strong. 



The major sources of uncertainty in our calculations 
is most likely to be connected with the treatment of 
the dust extinction and the assumed mass ra nge o f 
th e SN progenitors.^ In Dahlen fc Fransson ( | 1998 ) 



wei discuss xhese, and other possible sources or unccr 



tai hty. In this paper we also go deeper into the details 

of our model and include discussions on detectability, 
metal production and the influence of cosmological 
models on the derived rates. 



6. CONCLUSIONS 



We have found that the predicted number of core 
collapse SNe is rather insensitive to the assumed star 
formation scenario, as long as the star formation is 
calculated to match the same observed luminosity 
densities. Different star formation models may, how- 
ever, be tested by using high redshift subsamples of 
core collapse SNe. 

When it comes to the observable rates of Type la 
SNe, we find that these are more sensitive to the star 
formation and dust modeling. This is due to the 
fact that the Type la's are less linked to the envi- 
ronment where their progenitors were formed. The 
uncertainty in the life-time of the progenitors, com- 
bined with the sensitivity on the rates to the onset of 
star formation in monolithic- like models, contributes 
to the difficulty with using Type la's as probes for ei- 
ther different star formation scenarios, or progenitor 
models. 
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The number of core collapse SNe that could simulta- 
neously be detected with the NGST down to 1 nJy 
(mAB =31.4) in a 4x4 arcsec field is calculated to be 
~ 47 in the K' band. In order to detect SNe with 
z> 2, long wave-length filters (i.e. M') should be 
used. This is due to K-correction. The effect of the 
K-correction is enhanced by the spectral cutoff, caus- 
ing a majority of the SNe to drop out of the shorter 
wave-length filters at relatively lower redshifts. The 
high number predicted for the counts implies that 
searches for these SNe will be very important when 
it comes to the understanding of the star formation 
and the nucleosynthesis in the early Universe. 

The number of simultaneously detectable Type la 
SNe with the NGST limits are ^5-35, depending 
on model. Additional uncertainties widen this range 
even more. Of the simultaneously observable Type 
la's, only about 5% are at the rise of the light curve. 
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